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Abstract 

We study the phase diagram of S = 3/2 and S = 2 bond-alternating spin 
chains numerically. In previous papers, the phase diagram of S = 1 XXZ spin 
chain with bond-alternation was shown to reflect the hidden Z2 x Z2 symme- 
try. But for the higher S Heisenberg spin chain, the successive dimerization 
transition occurs, and for anisotropic spin chains the phase structure will be 
more colorful than the 5 = 1 case. Using recently developed methods, we 
show directly that the phase structure of the anisotropic spin chains relates 
to the Z2 x Z2 symmetry. 



1 



Typeset using REVTpX 



Haldane predicted that the antiferromagnetic Heisenberg spin chain has different behav- 
iors between integer and half odd integer spins [I]. For half odd integer spin cases, the ground 
state has gapless excitations and spin correlations decay in the power law. The ground state 
properties for half odd spin Heisenberg chains are the same ones of the 5=1/2 case. On 
the other hand for integer spin cases, the ground state is a singlet with finite excitation gap 
and the spin correlations decay exponentially. Another aspect for the Haldane gapped state 
was presented by Affleck et.al. [@]. They studied a special S = 1 SU(2) spin chain which 
has the property of the Haldane gap conjecture, and proposed valence bond solid (VBS) 
states. The S = 1 VBS state has the hidden antiferromagnetic structure. Later den Nijs 
and Rommelse found the string order parameter as the order parameter for it. 

Using a non-local unitary transformation, Kennedy and Tasaki [|J showed that the broken 
symmetry of S = 1 Haldane gap systems is a Z 2 x Z 2 type. Nearly four-fold degenerate 
states of the system with open boundary conditions reflect the discrete symmetry. For the 
higher S cases, Oshikawa || showed that odd integer S VBS states break the hidden Z 2 x Z 2 
symmetry, while even integer S VBS states do not violate the hidden Z 2 x Z 2 symmetry, 
and considered the generalized string order parameters. He pointed out another possibilities 
for the broken symmetry of higher S cases. 

In our previous papers with Okamoto HjTj, we studied the S = 1 XXZ spin chains with 
the bond-alternation, 

N 
3=1 

and we obtained the phase diagram reflecting the hidden Z 2 x Z 2 symmetry. Comparing 
with the phase structure of the quantum Ashkin- Teller model |8j which has a Z 2 x Z 2 
symmetry explicitly, the possible phase diagram consists of the two dimensional(2D) Ising, 
the Gaussian, and the Berezinskii-Kosterlitz-Thouless(BKT) transitions. Combining the 
quantum Ashkin- Teller picture with the generalized hidden Z 2 x Z 2 symmetry, we predicted 
that for -1 < A, -1 < S < 1, there are 2S + 1 2D Ising, 2S Gaussian, and 2S + 1 BKT 
critical lines. Especially the Gaussian lines should separate several VBS states, as was argued 
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by Guo et.al. and Oshikawa ||. According to them, the VBS state with the periodic 
boundary condition can be written by 

\S, M, PBC) 

= ^ N b[-bU\) S - M * 

N/2-1 

11 \ a 2j-l°2j °2j-l a 2j) \ a 2j°2j+l °2j a 2j+l) 

3=1 

x(a^ N -b^A) S+M \0), (2) 

where we describe the spin state by the Schwinger bosons, that is, creates the 5 = 1/2 

|(|) state at the j-th site, and M is an integer for integer S, or a half integer for half integer 
S, which changes from — S to S if we vary 5 from —1 to 1. Hereafter we denote this VBS 
state as M-VBS state. Oshikawa |J showed that for integer S case, the hidden Z 2 x Z 2 
symmetry is broken in the M-VBS state when S — M is an odd integer. For the isotropic case 
(A = 1), the 5 = 1 case is known that the transition between the Haldane and the dimer 
phases occurs about 5 = 0.26 pT0^P~5],[7[| . For the S = 3/2 case, Yajima and Takahashi |n 



studied with the density matrix renormalization group method and evaluated the transition 



point between M = 1/2 and M = 3/2 phase as #1/2-3/2 = 0.42 ± 0.02. Yamamoto fl3 
obatained consistent results as #1/2-3/2 = 0.43 ± 0.01 by quantum Monte Cairo calculation. 
For the S — 2 case, Yamanaka, Oshikawa, and Miyashita Jl5] evaluated the transition point 



by quantum Monte Carlo calculation. They obtained that the transition point between 
M = and M — 1 phases is in the region 0.05 < #0-1 < 0.3, and the transition point 



between M = 1 and M = 2 is in the region 0.5 < #i_2 < 0.6. Yamamoto [|T^] evaluated the 
transition point of S — 2 case more accurately as #0-1 = 0.18 ± 0.01, #i_2 = 0.545 ± 0.005. 

In this paper, we show the phase diagram of the model([L]) for S = 3/2 and 5 = 2 cases 
with the exact diagonalization method for finite size systems. In the following, we mainly 
argue the successive dimerization of the model(|T]). 



We use the following sine-Gordon Euclidean action as an effective model [IB] of the 
Hamiltonian(|l]) 
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where a is a lattice constat and t> is the sound velocity. The dual field 9 is defined as 



d 

VOT 

d_ 

dx 



■jL(iK9(T,x)), 



T, X 



d 
vdr 



(iK9(r,x)). 



We make the identification (J) = (f) + \^2ir, 9 = 9 J r \/2ix. For the free field theory, the primary 
operators exp(m\/2# + im\p2<$) have the scaling dimension x n ^ m = n 2 /2K + m 2 K/2 and 
the spin s n m = nm (where integer variables n and m are electric and magnetic charges in 
the Coulomb gas picture). The second term of eq.(f|) is the mass term for the Haldane gap 
systems |16,|17|]. (According to Affleck 0], we have yi oc cos 0/2, where is the topological 
angle.) 

After scaling a — > ae dl , we obtain the following 1-loop renormalization group equations 



d 1 
JlK 
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dyi f K\ 1 
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for K > 1. The third term 



From these equations we can neglect the third term of eq.( 
is important for the 2-D Ising transition and the level- 1 SU(2) WZW point, but not for 
the BKT and the Gaussian transitions. Neglecting y 2 we obtain the following Kosterlitz's 
recursion relation 

d 1 1 o 



dlK 8^' 



dyi f 



(5) 



Near K = 4, there occurs the BKT transition, and for 4 > K > 1 the Gaussian transition 
occurs at y± — 0. 

The obtained phase diagrams are in Fig.|I| for S = 3/2 case and in Fig.|2] for S = 2 case 
using the method in fT9| , |2!l| | . The topology of the phase diagrams is consistent with the 
prediction as is described above. 



At this stage, we will explain the method to determine the Gaussian line and the relation 
to the Z 2 x Z<i symmetry. Considering the hidden Z 2 x Z 2 symmetry, the successive dimerized 
transitions are of the Gaussian type. The Hamiltonian is invariant under the space inversion 
P:j —>■ N — 1 + j and the spin reversal T:a) <-> and all states(0) have the same eigenvalues 
of P and T (P = T = (— l) SN , independent of M). For the transition from M phase to 
the M + 1 one, the ground state is singlet in both phases, so we can not simply apply 
the phenomenological renormalization group technique. But with the twisted boundary 
condition 

the M-VBS state® changes to [0] 

\S, M, TBC) 

= (a] v 6 t 1 + 6] v a t 1 )^ M x 

JV/2-1 

11 \ a 2j-l u 2j u 2j-l a 2j) \ a 2j u 2j+l u 2j a 2j+l) 

3=1 

x(a^ N -b' N _A) S+M \0). (6) 

From this equation, the M-VBS state has P = T = (-i)Sn-s+m and the M + 1 state has 
P = T = (—\^sn-s+m+i^ so the two states have the different symmetries, and the energy 
levels cross at the transition point. 

In the sine-Gordon Language, the effect of the twisted boundary condition is to shift 



the magnetic charge m asm — > m + 1/2 [18|, so only the half integer magnetic charge 



operators exist. Then we have the following size dependence of the scaling dimensions of 
the m = ±1/2 operators a/2 cos(</>/\/2) [2^,1/2] and \/2 sin(0/\/2) [^0,1/2] EH' 



°' 1/2V J 2rrv 8 2 
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where L is the system size, y\ is a renormalized function of L as 



2n\ K / 2 - 2 



E e L q is the ground state energy of the periodic boundary conditions, and we denote E^° n 
as the corresponding energy to \[2 cos(0/v2), \/2sin(0/\/2) (e,o mean the parity P = 1, 
—1). In this method, the correction from the third term of eq.(||) does not affect since this 
method uses the duality of the quantum Ashkin- Teller model JF]. Note that Kolb |21| studied 



the XXZ spin chain with the twisted boundary condition S^ +1 = e ±% ®S^, S%- + i = S* and 
showed that for the half integer spin model the ground states are two-fold degenerate at 
$ = 7r, since for the half integer S XXZ spin chains, the second term in eq.(|J) is absent 
= 0). 

For S = 3/2 case, we study the N = 6,8, 10, 12 systems. In Fig. |, we show the two 
low lying energies of the subspace YJ S z = with the twisted boundary condition for A = 1, 
N = 12. For small S(> 0) the lowest state has P = T = 1, and for 1 > 5 > #1/2-3/2 the 
lowest state has P = T = — 1, as expected from eq.(^jj). Fig. [| shows the size dependence 
of the crossing point. We extrapolated the crossing point as polynomials of 1/N 2 |7| and 
the evaluated value is #1/2-3/2 = 0.4315, which is consistent with the previously obtained 



value ]1^,|T3[. For the BKT transitions between the XY and the VBS phases are determined 
by the similar method in ref |20| (see also 0,0]). The evaluated multi-critical points of 
the XY, M-VBS, and M + 1 VBS phases are (A M ,5 M ) = (0.8188,0), (0.5766, ±0.4052). 
For the transitions between the Neel and the VBS phases (2D Ising type), we use the 
phenomenological renormalization group method. 

For S = 2 case, we calculate the energy eigenvalues for N = 4,6, 8, 10. In Fig. ||, we show 
the two low lying energies of the subspace S z = with the twisted boundary condition 
for A = 1, N — 10. The level crossing occurs two times for < 5 < 1. For < 5 < #o_i 
the lowest state has P = T=1(5' = 2 Haldane phase, M — 0), for 5 -i < 5 < #i_2 the 
lowest state has P = T = — 1 (M = 1), and for 5\-2 < S < 1 the lowest state has P = T = 1 
(M = 2), as expected from eq.@. The hidden Z2 x Z2 symmetry is broken for M — ±1 
phase. Fig. || shows the size dependence of the crossing point (#0-1) • The extrapolated 



values of the transition points are 5o-i = 0.1830 (between M = and M = 1 VBS phases), 
and <5i_2 = 0.5505 (between M — 1 and M = 2 VBS phases). These values are consistent 



with the results of Yamamoto p3| . For 5 = 0, the transition point between the XY and the 
S = 2 Haldane phases is A c = 0.9650. The evaluated multi-critical points between the XY 
and the VBS phases are (A M ,5 M ) = (0.9708, ±0.1823) and (0.8401, ±0.5325). 

In this paper, we showed the phase diagrams of the S — 3/2, 2 bond-alternating XXZ spin 
chain. On the base of the hidden Zi x Z2 symmetry, we can see the successive dimerization 
numerically using the twisted boundary conditions. The two lowest excitations in the twisted 
boundary condition correspond to the order parameter and its dual disorder parameter (with 
respect to the Z 2 x Z 2 symmetry) |7|], so their crossing point is the self-dual point (in the 
Ashkin- Teller language) or the Gaussian critical point. Therefore, we can determine the 
critical point more precisely than the method using the string order parameter only ||12|| . 
especially near the BKT multicritical point. 
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FIG. 1. Phase diagram of the S = 3/2 model. The solid circles and the horizontal line separate 
M = -3/2,-1/2,1/2,3/2 VBS phases. Transitions between the XY and the VBS phases are of 
the BKT type. Transitions between the Neel and the VBS phases are of the 2D Ising type. 
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FIG. 2. Phase diagram of the S = 2 model in the region (a):-l < A < 1, and (b):0.8 < A < 1. 
Solid circles separate M = — 2,— 1,0, 1,2 VBS phases. 
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FIG. 4. Size dependence of the crossing points (#1/2-3/2) f° r S = 3/2 and A = 1. 
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FIG. 6. Size dependence of the crossing points (5q-i) for S = 2 and A = 1. 
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